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Abstract 


Highly  effective  treatments  for  the  effect  of  organophosphonate  and 
organophosphate  (OP)  nerve  agents  have  eluded  the  medical  community.  It  is 
known  that  acetylcholine  overload  is  one  of  the  effects  of  OP  toxicity,  but  the 
cellular  processes  leading  to  cardiac  toxicity  are  still  incompletely  imderstood. 
This  study  details  high  performance  computer  simulations  of  the 
electrophysiology  in  atrial  toxicity.  It  shows  that  hyperkalemia  of  the  tissue,  one 
of  the  manifestations  of  OP  intoxication,  promotes  the  processes  leading  to 
reentry,  a  precursor  of  atrial  fibrillation.  Then,  we  demonstrate  that  changes  in 
two  of  the  potassium  membrane  currents,  iicr  and  iKs  /  can  modulate  the  reentry 
process.  This  suggests  that  Class  III  anti-arrhythmic  agents  that  primarily  block 
these  currents  in  the  cardiac  cells  are  important  candidates  for  therapeutics  of  OP 
poisoning. 
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1.  Introduction 


Organophosphates  (OP)  and  organophosphonates'  constituents  of  insecticides 
and  nerve  agents,  respectively,  are  highly  toxic  substances  with  severe 
cardiovascular  ramification  upon  exposure.  An  estimated  95%  of  the  contacting 
agent  winds  up  within  and  is  transported  by  the  vasculature  and  is  primarily 
deposited  in  tissues  of  the  cardium.  For  a  review  of  its  effects,  see  Baskin  and 
Whitmer  [1],  Roth  et  al.  [2],  and  Hassler  et  al.  [3]. 

Toxicity  is  expressed  when  OP  binds  to  acetylcholinesterase  (AChE),  preventing 
hydrolysis  and  resulting  in  an  accumulation  of  acetylcholine  (ACh)  in  the  tissue. 
This  overload  causes  negative  inotropy  and  arrhythmia  not  only  in  the  atria 
(bradyarrhythmia),  but  also  in  the  Purkinje  fibers,  and  produces  electrical 
conduction  anomalies.  In  addition,  the  absence  of  cholinergic  agonists  increases 
the  calcitun  influx  and  stimulates  its  uptake  in  the  sarcoplasmic  reticulum  (SR), 
thereby  increasing  the  intracellular  calcium.  Accompanying  this  process  is  an 
increase  in  the  vagal  tone  that  may  contribute  to  tiie  generation  of  atrial 
fibrillation.  Chohnergic  stimulation  of  atrial  tissue  reproducibly  induces  atrial 
fibriHation  (AF)  [4-7]. 

The  physiological  events  that  culminate  in  AF  upon  cholinergic  stimulation  are 
only  incompletely  imderstood.  Some  of  the  known  factors  ancillary  to  AF  are 

(1)  non-uniform  shortening  of  the  myocardial  refractory  period  [7], 

(2)  production  of  intra-atrial  conduction  delays, 

(3)  bradycardia  that  increases  susceptibility  to  reentry, 

(4)  shortening  of  the  atrial  action  potential  duration, 

(5)  cholinergic  stimulation  may  be  a  sine  non  quo  for  the  preparation  of  the 
subtrate  where  atrial  reentry  may  occur, 

(6)  inhomogeneous  repolarization  may  be  necessary  for  atrial  arrhythmia,  and 

(7)  heterogeneity  of  the  substrate. 

OP  toxicity  is  accompanied  by  elevation  of  the  concentration  of  [K+]o,  the 
extracellular  potassium  concentration  that  influences  the  transmembrane 
resistance  in  Purkinje  fibers,  and  [Na*]o,  whose  elevation  mimics  ischaemia.  In 
the  atria,  a  biphasic  behavior  of  conduction  velocity  as  a  fvmction  of  potassium 
concentration  has  been  noted.  Moderate  elevation  of  extracellular  potassium  (up 
to  8  mM)  increases  the  conduction  velocity;  but  in  the  concentration  range  of 
8-16  mM,  slowdown  and  eventual  termination  of  conduction  has  been  observed. 
Changes  in  the  action  potential  due  to  toxicology-caused  shifts  in  potassium 
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concentration  also  form  the  basis  of  an  updated  model  of  soman-caused  cardiac 
toxicity  [8]. 

The  tremendous  increase  in  computational  resources,  especially  parallel 
processing  coupled  with  novel  experimental  techniques,  has  allowed  the 
emergence  of  models  of  subceUular  dynamics.  Thus,  it  is  possible  to  emulate 
cellular  processes  and  study  the  effect  of  the  presence  of  toxic  agents  at  a 
fundamental  level  effectively  and  without  animal  experimentation. 

The  following  sections  describe  a  computer  study  of  the  effect  of  ACh  overload 
on  the  electrophysiology  of  atrial  tissue,  especially  the  processes  leading  to 
reentry,  that  subsequently  and  usually  culminate  in  AF,  the  precursor  of  sudden 
cardiac  death  (SCD). 


2.  Materials  and  Methods 


This  report  details  our  approach  to  the  simulation  of  cardiac  toxicity  created  by 
the  presence  of  OP  in  atrial  tissue.  This  study  especially  seeks  to  identify 
conditions  that  promote  the  modulation  of  membrane  currents  that  may  lead  to 
reentry,  a  precursor  of  AF  and  ventricular  fibrillation  (VF).  Conversely,  it  seeks 
to  answer  the  question  of  whether  it  is  possible  to  identify  conditions  whose 
modulation  may  reverse  the  process  and  thus  suggest  strategies  for  therapies. 
The  approach  was  to  simulate  conditions  that  result  from  OP  intoxication  and  to 
see  if  these  conditions  favor  reentry  of  the  action  potential. 

2.1  Prior  Work 

Two  models  of  AF  have  been  advanced  — the  rapidly  discharging  focus  and 
reentrant  excitation.  Here  we  discuss  the  latter,  where  atrial  fibrillation  is 
manifested  by  rapid,  irregular  action  potential  waves  of  variable  wavelength  and 
timing  [9-11]  along  reentrant  paths  extinguishing  that  are  also  being 
reestablished.  These  irregular  waves  prevent  the  proper  functioning  of  the  atria 
due  to  the  shortening  of  the  refractory  period  and  the  loss  of  the  normal 
lengthening  of  refractoriness  when  the  heart  slows  down.  The  wave  of  excitation 
returns  to  die  tissue  it  previously  exited  and  reexcites  it,  hence  the  name  reentry. 
Until  quite  recently,  it  was  believed  that  below  a  certain  volume  of  atrial  tissue 
reentry,  one  of  the  forms  of  AF  could  not  be  sustained.  However,  the  work  of 
Wakimoto  et  al.  [12]  showed  that  even  in  mtnine  atria,  AF  is  inducible.  The 
breakup  of  the  spiral  waves  signals  the  transition  from  tachycardia  to  fibrillation. 

Moe  and  his  collaborators  [13]  were  the  first  to  give  a  mathematical  description 
of  atrial  fibrillation.  Several  electrophysiological  models  of  the  atria,  such  as 
those  of  Ramirez  et  al.  [14],  Lindblad  et  al.  [15],  and  Nygren  et  al.  [16],  based  on 
the  Hodgkin  and  Huxley  formahsm  and  patch  clamp  data,  are  now  available. 
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These  models  allow  in  silico*  studies  of  the  effect  of  changes  in  the  membrane 
currents  caused  by  pathological  conditions  and  their  effect  on  the  subsequent 
changes  in  the  action  potential  to  be  studied  without  resort  to  in  vivo 
experiments. 

Table  1,  from  Koumi  et  al.  [17],  gives  an  indication  of  the  effect  of  ACh  on  atrial 
myocytes.  Significant  (38%)  decrease  in  the  action  potential  duration  was  noted 
when  1  pM  of  ACh  was  applied  to  the  atrial  tissue.  ACh  also  increases 
conductance  of  the  tissue  and  especially  that  of  the  Ikt  channel  which  occurs 
during  myocardial  ischemia. 

Table  1.  Effect  of  ACh  on  atrial  myocytes  from  Koumi  et  al.  [17]. 


Resting  Membrane 
Potential 
(mV) 

Action  Potential 
Amplitude 
(mV) 

Action  Potential 

Duration  at  90% 
Repolarization 
(ms) 

Control 

-73.0  ±  7.2 

16.8  ±  5.0 

289.4  ±23.0 

1  pMACh 

-77.4  +  6.1 

-8.1  ±  4.8 

112.6  ±  14.5 

2.2  The  Simulation  Setup 

For  these  calculations,  an  atrial  tissue,  3x3  cm,  and  ubsequently  one  of 
3x4  cm,  with  its  electrophysiology  described  by  the  Nygren  et  al.  [16]  model, 
was  stimulated  by  a  pulse  of  40  pA/cm^  of  2  ms  in  duration  at  one  edge.  The 
presence  of  OPs  was  modeled  by  patches  of  tissue  having  a  [K+]o  of  10.8  up  to 
25.0  mM  instead  of  the  ambient  5.4  mM.  This  expressed  the  observation  that 
h5q)erkalemia  is  one  manifestation  of  OP  toxicity.  In  these  regions,  the  resting 
membrane  potential  is  less  negative,  and  depends  on  the  potassium 
concentration.  Admittedly,  this  is  coarse  graining  the  problem  since  localized 
sodium  and  calcium  concentration  deviations,  among  other  changes  found  in  the 
presence  of  OP,  are  also  indicators  of  the  change  of  the  state  of  the  cell.  For  these 
calculations,  the  tissue  was  represented  as  a  slab  of  300  x  300  x  1  nodes  and  also 
as  300  X  400  x  1  with  node  spacing  of  0.01  cm.  The  pacing  used  to  establish  the 
initial  conditions  consisted  of  ten  10  pulses  at  300-ms  intervals. 

2.3  Computational  Approach 

The  computational  experiments  reported  were  designed  to  establish  and 
simulate  conditions  that  favor  the  establishment  of  reentry  in  atrial  tissue 
affected  by  OP. 


*  In  silico  refers  to  the  modeling  of  living  processes  on  the  computer. 
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Cardiac  tissue  is  usually  modeled  as  a  bidomain,  as  two  interpenetrating 
domains  representing  the  inside  and  the  outside  of  the  cell,  separated  by  a 
membrane  through  which  the  current  transits  from  one  domain  to  the  other.  In 
each  of  these  domains,  the  potential  is  described  by 


(1) 


and 


=  (2) 

Here,  the  <p,  and  ^  are  intracellular  and  interstitial  potentials,  D's  are  the 
diffusion  tensors,  and  the  ^  is  the  cellular  surface  to  volume  ratio.  Im  is  the 
membrane  current,  and  hi  and  he  are  the  current  source  densities.  When  the 
diffusion  tensors  of  the  two  domains  are  proportional  to  each  other,  as  was 
assumed  here,  the  bidomain  equations  simplify  to  the  monodomain  model. 
Monodomain  calculations  are  considerably  less  expensive  to  perform  than 
bidomain  calculations.  For  the  calculations  reported  here,  fiber  orientation  (one 
of  the  diffusion  matrix  entries)  was  assumed  to  be  uniform. 

Two  other  sets  of  equations  need  to  be  solved:  equations  describing  the  gating  of 
the  ion  channels  and  equations  for  the  concentration  changes  of  the  ions. 
Generically,  these  may  be  written  as  follows: 

dyidt  =  (y„  -  ,  (3) 

where  y  represents  a  gating  variable  and  Zy  the  time  constant.  In  addition,  the 
ion  concentrations  had  to  be  followed.  For  calcium,  for  example,  the  equation  is 
as  follows: 

d[Ca\ldt  =  +0.07(10"^  -[Cfl],).  (4) 

In  equation  (4),  the  quantities  in  the  square  brackets  are  the  ion  concentrations 
and  hi  represents  the  L-type  Ca++  current.  Finally,  the  individual  membrane 
currents  from  experimental  measurements  are  included.  No-flux  boundary 
conditions  were  enforced  at  the  edge  of  the  tissue. 

The  solution  consists  of  three  steps:  (1)  calculation  of  the  ionic  currents,  Im,  (2)  the 
determination  of  the  transmembrane  potential  Vm,  the  difference  between  the 
intra  and  the  interstitial  potentials,  and  (3)  the  upgrading  of  the  field  potentials. 
In  step  (3)  conjugate-gradient  was  the  solution  method  used. 

Equal  grid  size,  0.01  cm  in  both  directions,  was  used.  Some  experimentation 
with  time  step  size  was  required,  but  0.001  ms  yielded  stable  results. 
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The  calculations  were  carried  out  with  the  code  CardioWave  on  the  computer 
assets  of  the  Major  Shared  Resource  Center  (MSRC)  at  Aberdeen  Proving 
Gromid,  MD. 

In  parallel  mode,  a  typical  rrm  on  16  nodes  on  a  SGI  Origin  3000  required  in 
excess  of  7  hr  (in  non-dedicated  mode)  for  a  grid  of  300  x  300  to  simulate  300-ms 
real  time.  For  a  larger  number  of  processors  improvement  in  wall  clock  time  was 
achieved,  but  the  speedup  was  not  absolutely  linear,  as  seen  in  Figure  1. 


Figure  1.  Performance  of  the  code  using  the  shared-memory  Origin  3000. 


3.  Results 


We  report  on  three  computational  experiments. 

(1)  A  pulse  of  40  |iA/cm2  of  2  ms  duration,  starting  at  the  left  edge,  propagates 
into  the  tissue,  where  it  encovmters  two  subregions  of  higher  than  normal 
potassium  concentration.  Each  of  these  subregions  was  1  cm2  size.  The 
first  was  located  along  the  central  axis,  flush  with  the  left-hand  bomidary. 
The  second,  of  the  same  size,  was  located  along  the  lower  edge  of  the 
tissue,  1  cm  from  the  left  boundary. 
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(2)  The  second  in  silica  experiment  consisted  of  the  same  tissue  with  one  region 
of  higher  potassium  concentration  along  the  central  axis  of  the  tissue,  1  cm^ 
in  size,  located  at  the  left-hand  boundary,  aligned  normally  to  the  action 
potential  advancing  from  the  left-hand  boundary. 

(3)  Finally,  the  same  field  of  the  action  potential  used  as  in  (2)  was  examined, 
but  with  the  iio-  and  ixs  currents  blocked  to  99%  of  its  nominal  value.  This 
was  based  on  the  insight  that  new  Class  111  antiarrhythmic  drugs,  such  as 
Azimilide,  are  potassium  channel  antagonists.  Modulating  these  currents 
prolongs  the  duration  of  the  action  potential  and  the  effective  refractory 
period.  It  was  felt  that  this  antiarrhythmic  characteristic  may  translate  into 
modulating  the  processes  that  favor  reentry  since  it  is  known  that  this  drug 
reduces  the  incidence  of  AF  and  atrial  flutter. 

Figure  2,  panels  (a-d)  show  the  encounter  of  the  action  potential,  started  at  the 
left  edge  of  the  atrial  tissue  as  it  encounters  a  subregion  with  a  potassium 
concentration,  double  the  nominal  value.  This  mimics  cardiac  toxicity  caused  by 
the  presence  of  an  OP-based  substance. 

In  Figure  3,  the  action  potential  speeds  up  in  the  region  with  the  higher  than 
normal  potassium  concentration.  The  curvature  m  the  action  potential  wave 
front  is  accentuated  as  time  progresses.  Behind  and  directly  in  line  with  the 
concentration,  Vm,  the  voltage,  continues  to  drop  sharply  in  contrast  to  the 
adjacent  areas  where  the  drop  is  more  gradual  and  in  line  with  healthy  tissue 
dynamics. 

The  front  continues  to  distort,  and  aft  of  the  heterogeneous  region  boundary 
undulations  in  the  voltage  equipotential  lines  appear.  The  front  is  no  longer 
coherent,  and  the  breakup  of  the  integrity  of  the  wave  has  commenced. 

The  second  set  of  computer  experiments  with  multiple  regions  of  higher  than 
normal  potassium  concentrations  shows  similar  evolution  of  the  wave.  The  time 
history  of  the  action  potential  at  a  particular  location  next  to  the  high  potassium 
concentration  mirrors  the  trend  toward  disorgaruzation  in  the  wave  structure, 
signaling  eventual  breakup. 

In  the  final  Figure  4,  the  iKi  and  ixs  currents  have  been  reduced,  mimicking  the 
effect  of  antiarrhythmic  drug  treatment  for  AF.  Lengthening  the  action  potential 
duration  (APD)  is  one  of  its  effects.  This  effect  is  also  apparent  when  contrasting 
Figures  3  and  4.  The  change  in  behavior  is  a  clear  indicator  that  the  search  for  a 
therapy  for  OP  toxicity  needs  to  be  investigated  because  pharmaceutics  that 
modulate  the  potassium  currents  prolong  APD,  thus  preventing  reentry. 
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(c)  (d) 

Figure  2.  Frame  sequence  showing  the  evolution  of  the  action  potential  in  atrial  tissue  using  the 
Nygren  model  for  the  electrophysiology  of  the  tissue.  The  times  of  the  frames  going 
clockwise  are  (a)  25  ms,  (b)  60  ms,  (c)  140  ms,  and  (d)  200  ms.  The  stimulus  consisted 
of  a  current  pulse  of  40  |aA/  cm^  delivered  to  the  left-most  column  of  atrial  cells.  The 
stimulus  lasted  for  2  ms;  VM  in  the  legend  is  the  voltage  in  mV. 
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Figure  3.  The  action  potential  in  an  atrial  tissue  containing  a  single  hyperkalemic, 
cardioplegic  subregion. 
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Figure  4.  The  effect  of  blocking  iKr  and  iKs  currents  on  the  action  potential. 


4.  Discussion 


In  this  study  we  examined  the  conditions,  in  terms  of  membrane  currents,  that 
promote  atrial  reentry,  the  precursor  of  the  processes  leading  to  AF.  One  of  the 
underlying  assumptions  was  that  OP-induced  ACh  overload  results  in  patches  of 
tissue  with  [K+]o  in  excess  of  normal  values.  When  an  action  potential  encounters 
such  a  region,  processes  leading  to  reentry  are  commenced.  Conditions  through 
the  modulation  of  iK  currents  were  sought  that  could  antagonize  reentry.  This  is 
important,  since  it  is  believed  that  the  morphology  of  AF  is  a  progression  from 
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tachycardia,  through  the  disintegration  of  spiral  waves  in  the  atrial  tissue,  to 
fibrillation. 

4.1  OP  Toxicity  Modeled  as  Extracellular  Potassium  Concentration 

Changes  in  ion  movement  are  central  to  the  expression  of  OP-induced  cardiac 
toxicity.  An  ancillary  aspect  is  that  the  accumulation  of  ions  in  the  interstitial 
space  may  play  an  important  role  in  the  overall  behavior.  The  reasoning  is  as 
follows:  excessive  parasympathetic  stimulation  results  in  excessive  ACh  release. 
ACh  binds  to  M2  muscarinic  receptors,  activates  the  G-protein,  and  contributes  to 
the  hyperpolarization  of  the  the  SA  node.  When  ACh  activates  iicACh/  K*  exits  the 
cell.  The  OP-induced  ACh  overload  results  in  more  receptors  than  usual  being 
boimd,  leading  also  to  an  increase  in  the  external  concentration  of  potassium 
ions. 

At  high  [K*]o  in  the  tissue,  similar  to  the  condition  of  ischemia  or  during 
strenuous  exercise,  opening  less  than  1%  of  Ik^atp  channels  causes  significant 
shortening  of  AP. 

Excessive  concentration  of  potassium  ions  changes  the  concentration  gradients 
and  the  response  of  the  Na-K  pump.  In  the  initial  stages  of  the  accumulation  of 
the  potassium  ions,  the  cells  are  more  excitable,  but  as  the  membrane  potential 
drops  below  the  excitation  threshold,  the  generation  of  an  action  potential  is  no 
longer  possible.  Hyperkalemia  suppresses  automaticity  of  the  heart,  slows 
conduction,  and  leads  to  bradycardia. 

4.2  Reentry 

It  was  shown  that  the  process  that  can  culminate  in  reentry  in  atrial  tissue  when 
a  supra-threshold  stimulus  40  pA/cm^  was  appHed,  since  propagates  into  an 
OP-modified  heterogeneous  atrial  tissue.  The  wave  front  that  resulted  was  a 
spiral  due  to  the  refractory  state  of  the  tissue.  Reentry  is  characterized  by  the 
action  potential  that  instead  of  transiting  and  exiting  the  tissue,  reenters 
previously  traversed  tissue. 

With  an  ectopic  beat  present  in  the  tissue,  the  encotmter  between  the  original 
stimulus  and  the  stimulas  representing  the  ectopic  event  produces  a  turning 
wave  front,  a  precursor  for  reentry. 

4.3  Modulation  of  Reentry  in  the  Atria 

Qu  et  al.  [18]  noted  in  their  study  of  reentry  in  ventricular  tissue  that  the 
condition  for  spiral  wave  breakup,  and  thus  the  condition  for  VF,  was  reached 
when  the  slope  of  the  action  potential  restitution  curve  exceeded  one.  Flattening 
the  restitution  curve  can  exert  an  anti-fibiillatory  effect. 
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In  the  atrial  tissue  studies  reported  here,  iicr  and  iKs  currents  were  modulated.  At 
1%  of  nominal  value,  considerable  modulation  of  the  action  potential  d}mamics 
is  observed.  The  block  of  the  inward  current  can  contribute  to  APD 
prolongation,  especially  at  shorter  cycle  lengths.  This  suggests  that  Class  III 
drugs  that  prolong  APD  have  the  potential  for  prophylaxis  or  treatment  of 
OP-caused  atrial  arrhythmias. 

4.4  Summary 

Computer  experiments  to  simulate  reentry  in  atrial  tissue  with  OP-induced  ACh 
overload  were  performed.  The  presence  of  the  OP  was  modeled  as  local 
potassium  overload.  The  experiments  showed  the  incipient  processes  that  can 
enhance  transition  from  atrial  tachycardia  to  atrial  fibrillation.  It  was 
demonstrated  that  modulation  of  membrane  current  affects  changes  in  the 
process  leading  to  reentry.  Possible  conditions  for  the  elimination  for  reentry  in 
this  tissue  model  were  found.  This  suggests  viable  new  approaches  to  preventing 
and  treating  OP  toxicity.  Further  computer  studies  and  planned  in  vivo 
experiments  at  the  U.S.  Army  Medical  Research  Institute  of  Chemical  Defense 
(USAMRICD),  Aberdeen  Proving  Grormd,  MD,  in  the  near  future  should  help  to 
validate  this  approach. 
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